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Abstract 

We investigate the propagation of a massless scalar field on a star graph, 
modeling the junction of n quantum wires. The vertex of the graph is rep- 
resented by a point-like impurity (defect), characterized by a one-body scat- 
tering matrix. The general case of off-critical scattering matrix with bound 
and/or antibound states is considered. We demonstrate that the contribution 
of these states to the scalar field is fixed by causality (local commutativity) , 
which is the key point of our investigation. Two different regimes of the theory 
emerge at this stage. If bound sates are absent, the energy is conserved and the 
theory admits unitary time evolution. The behavior changes if bound states 
are present, because each such state generates a kind of damped harmonic os- 
cillator in the spectrum of the field. These oscillators lead to the breakdown 
of time translation invariance. We study in both regimes the electromagnetic 
conductance of the Luttinger liquid on the quantum wire junction. We derive 
an explicit expression for the conductance in terms of the scattering matrix 
and show that antibound and bound states have a different impact, giving 
raise to oscillations with exponentially damped and growing amplitudes re- 
spectively. 
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1 Introduction 



Quantum field theory on graphs attracts recently much attention. Besides the purely 
theoretical interest, this attention is largely motivated by the physics of quantum 
wires. In fact, in many physical applications a quantum wire can be fairly well 
approximated by a graph. Each vertex of such a graph represents a junction of 
some wires. From the theoretical point of view the junction can be interpreted as a 
point-like defect, which is characterized by a scattering matrix S. The goal of this 
paper is to explore the case when S admits bound and/or antibound states. It is 
worth stressing that such states are absent at the scale-invariant (critical) points. In 
order to keep their contribution, one must explore therefore the theory away from 
criticality. In this respect the present investigation extends some previous studies 
P]-[T5] in the subject, which have been mainly focused on the scale-invariant limit 
of the theory. 

We consider in this paper graphs T of the type shown in Fig.l. They are called 
star graphs and represent the building blocks for generic graphs. We start by estab- 
lishing the analytic structure of the scattering matrix S associated with the vertex of 
T. In general S admits both bound and antibound states. We demonstrate that their 
contribution to the massless scalar field <p is fixed by causality (local commutativity) 
up to a free parameter t m . Two different regimes of the theory emerge at this stage. 
In absence of bound states, the field cp admits a unitary time evolution respecting 
time-translation invariance. The situation radically changes when S admits bound 
states. Each of them generates in the spectrum of the theory a kind of damped 
harmonic oscillator. These oscillators lead to a breakdown of time-translation in- 
variance. The energy of the system is no longer conserved, which signals a nontrivial 
flow of energy trough the boundary, being in this case the vertex of T. For t < t m the 
vacuum energy flow is outgoing and the system is in a dissipative regime. For t > t m 
the energy flow is incoming and the vacuum energy is growing. Although Hermitian, 
the underlying Hamiltonian is not self-adjoint, implying nonunitary time evolution. 
In what follows we refer to the bound states of S as Boundary Bound States (BBS), 
because they decay exponentially in the bulk of the graph. We concentrate on the 
massless scalar field <p because it is the fundamental block of bosonization. 

In the second part of the paper we extend the bosonization procedure from the 
line to star graphs and analyze the impact of the off-critical boundary conditions on 
the Luttinger liquid. The renewed interest in Luttinger liquids arises from the recent 
progress in producing narrow quantum wires with few, or even single, conducting 
channels. Examples are quantum Hall edges [IB], polydiacetylene [T7] and carbon 
nanotubes [TBI [EH] • As a potential application of our framework, we investigate the 
electromagnetic conductance of the Luttinger liquid on T away from criticality. We 



1 




Figure 1: A star graph T with n edges. 

derive an explicit expression of the conductance in terms of the scattering matrix 
and discuss in detail the different impact of bound and antibound states. 

The paper is organized as follows. In section 2 we describe the two regimes of 
a massless scalar ip field on T with off-critical boundary conditions at the vertex. 
We discuss both the scattering and BBS contributions. Special attention is paid 
to local commutativity. We introduce here also the dual field tp and the vertex op- 
erators defined in terms of ip and ip. In section 3 we derive the basic correlation 
functions. The vacuum instability due to the BBS and the accompanying break- 
down of time-translation invariance are also analyzed there. In Section 4 we study 
the Luttinger liquid on F and solve the Tomonaga-Luttinger model away from crit- 
icality. Coupling this model to a time-dependent electromagnetic field, we derive 
the conductance tensor in presence of (anti)bound states in S. We also show here 
the bound and antibound states imply different time behavior for the conductance, 
which represents a potential experimental signature for studying the analytic struc- 
ture of the scattering matrix at the vertex of T. Section 5 collects our conclusions. 
In the appendices we give some computational details and prove some statements, 
used in the body of the paper. 

2 Massless bosons with off-critical boundary con- 
ditions on a star graph 

2.1 Preliminaries 

We consider in this section the equation 

(d t -d 2 x )<p(t,x,i) = 0, x>0, i = l,...,n. (2.1) 

where t is the time and (x > 0,i = 1, ...,n) are the coordinates of any point P on 
the star graph T, shown in Fig. [TJ In order to solve (I2.ip . we must fix the boundary 
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condition in the vertex (junction) V of T. A standard QFT assumption at this point 
is to require that the operator K = — d 2 x on T is self-adjoint^. This requirement, 
which is trivially satisfied on the line R, is delicate on T. Fortunately however, 
the problem has been intensively investigated in the recent mathematical literature 
[2D]-[24J, where the subject goes under the name of "quantum graphs" [25]. From 
these studies one infers that K is self-adjoint on V if and only if the field (p satisfies 
the boundary condition 

n 

[A(I - U)ij <p(t, 0, j) - i(I + U) t ,(d x p)(t, 0,j)} = , Vt e R, i = 1, ...,n , 

3=1 

(2.2) 

where U is any unitary matrix and A > is a parameter with dimension of mass 
introduced in order to recover the correct physical dimensions. 

Eq. (12. 2p generalizes to the graph F the mixed (Robin) boundary condition on 
the half-line R + . The matrices U = I and U = —I define the Neumann and Dirichlet 
boundary conditions respectively. The physical interpretation of (I2.2p in the context 
of bosonization is established in [8] . At criticality ( 12. 2 j) describes the splitting of the 
[/(l)-charge current at the junction x = 0. We show below (see eqs. ( I4.20p . ( I4.22p ) 
that (12. 2p can be reformulated in terms of currents also away from criticality. 

Besides the boundary conditions, we must fix also the initial conditions imposing 
the equal-time commutation relations 

[cp(0, x u in) , <p(0, x 2 , i 2 )\ = [(<9^)(0, xi, i x ) , (d t <p)(0, x 2 , i 2 )} = , (2.3) 

[{d t (f)(0, xi.ii) , 95(0, x 2 , i 2 )\ = -iS ilia S(xi - x 2 ) . (2.4) 

As already mentioned in the introduction, for the explicit construction of (p it is 
convenient to interpret [26J the vertex of T as a point-like impurity (defect) [27]- [29], 
characterized by a nontrivial scattering matrix S. The S'-matrix is associated [20J 
to the operator K on Y and is fully determined by the boundary conditions (12. 2p . 
The explicit form of S is [20] [22] 

S(k) = -[A(I -U) + k(I + U)]- 1 ^! -U)- k(I + U)\ (2.5) 

and has a transparent physical meaning: the diagonal element Su(k) represents the 
reflection amplitude from the defect on the edge whereas Sij(k) with i 7^ j equals 
the transmission amplitude from E t to Ej. 
By construction (I2.5P is unitary 

S (k)* = S{k)- 1 , (2.6) 

5 In order to keep the BBS contribution we do not impose K > 0. 
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and satisfies Hermitian analyticity 

S(k)* = S{-k) . (2.7) 
Combining ( 12. 6 p and ( 12. 7p . one gets 

S(k) S(—k) — I. (2.8) 

Notice also that 

S(X) = U. (2.9) 

We see that the boundary condition (12. 2p is fixed actually by the scale A and the 
value of scattering matrix at that scale. 

Let us establish now the analyticity properties of (12. 5p in the complex fc-plane. 
Let U be the unitary matrix diagonalizing U and let us parametrize 

U d = U- l UU (2.10) 

as follows 

U d = diag (e 2iai ,e 2iQ2 ,...,e 2iQn ) , a; el. (2.11) 
Using (I2.5p . one easily verifies that U diagonalizes also S{k) for any k and that 

Si( k)=U^S(k)U = ^{^^,...,\±^] , (2.12) 

\k — 177! k — 1?72 k-i?7 n / 

where 

7T 7T 

rji = Atan(aj) , --<«»<-• (2.13) 

We have therefore a simple direct proof of the following statement [23] : 

Proposition 1: The scattering matrix S(k) given by $2.5\) is a meromorphic 
function with poles located on the imaginary axis and different from 0. 

The poles belonging to the lower and upper half-plane are known [30] as anti- 
bound and bound states respectively. It is natural to represent at this stage the field 
ip as a linear combination 

(p(t, x, i) = ip( s \t, x, i) + ip {b) (t, x, i) , (2.14) 

where collects the contribution of the scattering states and (f^ b \ that of its bound 
states. 
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2.2 The scattering component (p^ 

The scattering component <p s ) is known from previous studies [26]. One has 

<pM(t,x,i) = [°° d ^_ [a?(fc)e i (' fc l t - te ) + a i (fc)e- i (l fc l*- te )] , (2.15) 
i-oo 27Ta/2|A;| 



where {ai(fc), a*(fc) : k G M} generate the reflection-transmission (boundary) al- 
gebra [291 El] corresponding to the scattering matrix f |2.5j) . This is an associative 
algebra ^4 with identity element 1, whose generators {cfj(&), a* l {k) : k G M} satisfy 
the commutation relations 

a^fci) a i2 (fc 2 ) - ai 2 (fc 2 ) «u(^i) = , (2.16) 
a^fcx) a* 2 (k 2 ) - a? (fc 2 ) o? (fcx) = , (2.17) 
a h (ki) a* 2 (k 2 ) - a* 2 (k 2 ) a^fa) = 2tt [5 ili2 S(kx - k 2 ) + S ili2 (ki)S(k 1 + k 2 )} 1 , 

(2.18) 

and the constraints^) 

n n 

Oi(k) = ^S ij (k)a j (-k) , a K k ) = J2 a *j(- k ) S Ji(- k ) • ( 2 - 19 ) 

3=1 3=1 

Time reversal is implemented in the algebra A by 

Ta(k)T- 1 = a{-k) , Ta*{k)T~ l = a*(-k) , (2.20) 

T being antiunitary. This action is consistent with f l2.16tl2TT9|) provided that 

S{k) = S{-k), (2.21) 

where the bar stands for complex conjugation. Combining (12.71) and (12.211) . one 
infers^ 

S{kY = S{k) , (2.22) 

which, in view of (12. 9ft . implies 

U t = U. (2.23) 

Therefore, the boundary conditions (I2.2p respect time-reversal invariance only for 
symmetric U, which will be our choice in what follows. Accordingly, the diagonal- 
izing matrix U in (12.101) can be chosen to satisfy 

U 1 =W X , U = U. (2.24) 



6 Eq. (I2.8[) guarantees their consistency. 
7 The superscript t denotes transposition. 
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The case TP 7^ U with broken time reversal has been investigated recently in [IB] . 

The field (12.1 5p obviously satisfies the equation of motion ( 12. ip . Using (I2.19p . 
one can show after some algebra that (p^ satisfies the boundary condition ( 12. 2 p as 
well. By means of fl2.16H2.i8]) and fl 2 . 2 2 [) one can also verify that <p( s > satisfies the 
canonical relations f 12 . 3 j) . As far as ( 12 .4p is concerned, one obtains 

/°° Ah 
—e ik ^S tll2 (k), 
-00 ^ 

(2.25) 

with X12 = x\ + x 2 . The integral in the right hand side of ( I2.25P can be computed 
using the analytic properties of S(k). Let us denote by V+ = {irj : rj > 0} and by 
V- = {\rj : r\ < 0} the poles of S(k) belonging to the upper and lower half-plane 
respectively. V + collects the bound states, whereas P_ the antibound states. The 
set V + (V-) coincides with the distinct positive (negative) rji appearing in eq. (12. 13j) . 
Let us introduce also the matrix 

B^l = - lim (* - iv)S ili2 (k) , ir] G V± . (2.26) 

Now, keeping in mind that xu > 0, the integral in (I2.25P can be computed by means 
of the Cauchy integral formula. One finds 

[(V s) )(0,x 1 ,^ 1 ),^(0,x 2 ,^)] = -i^ 2 5(x 1 -a; 2 ) + i ^ rje'^ R%1 , (2.27) 

iri£V+ 

Let us briefly analyze this result. 

If V + = (no poles in the upper half-plane) the second term in the right hand 
side of (12.271) vanishes and the field tp^ satisfies the canonical commutation relations 
(12.31 I2.4p . Accordingly, we set in this case cp^(t,x,i) = 0. This is the regime 
previously investigated in [TIE]- 

If V + 7^ 0, the scattering matrix S(k) admits bound states. The scattering states 
(represented by the plane waves in (I2.15P ) do not form a complete set, which leads to 
the second term in the right hand side of (I2.27p . In order to compensate this term, 
one has to introduce new boundary degrees of freedom collected in (p(°)(t,x,i) 7^ 0. 
Since the guiding line for constructing ^ is local commutativity, it is instructive 
to compute the commutator of ip^ at generic points on the graph. One finds 

[ip {s \t!,x u ii) , <p {s) (t 2 ,x 2 ,i2)} = 

~ [e(*i2 + X&) + e{t 12 - x 12 )] 6 ilia - - [e(t 12 + x 12 ) + e(t 12 - x 12 )} S ili2 {0) 

~\ E + 5 ? i2)e-" {tl2+Sl2) - 0(-h 2 + x 12 )e^-^} R^i 
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~\ E \- 6{t ^ - ^)e r ' ( * 12 - Sl2) - 0{-t l2 - x 12 )e~^ + ^} R%1 , (2.28) 

where e is the sign function, t\ 2 = h — t 2 , %i2 = x± — x 2 and X12 = x± + x 2 . From 
( I2.28P one gets at space-like separations t\ 2 — x\ 2 < 

[tpUfaxuH) , ^ s \t 2 ,x 2l i 2 )} = i ^ e"^ 2 sinhfe) R^ 2 . (2.29) 

Summarizing, we proved the following statement. 

Proposition 2: The scattering component (p^ is a canonical local field only in 
absence of boundary bound states (V+ — 0). 

In the case V+ 7^ our strategy below will be construct a boundary component 
(p( b \ such that the total field if = ip^ + ip( b > is both canonical and local. 

2.3 The BBS component tp® 

We assume in this subsection V + 7^ and introduce the set 1+ = {i : rji > 0}, where 
r]i are given by ( 12. 13ft . The wave functions of the bound states of K = —d 2 x on F 
are then {e~ ViX : i 6 They are normalizable and generate the following set 

^ e -Vt(x±t) . ^ e j + j Q f rea j so i u tions of the equation of motion ( 12. ip . We associate 
with each index i G X + a quantum oscillator. The quantum boundary degrees of 
freedom are described therefore by an algebra B, generated by {h,b* : i G 
which satisfy 

&i> &i 2 - b l2 b tl =0, 6? b* 2 - b* 2 bl=0, (2.30) 
K K - K b h = 5 ni2 (2.31) 
and commute with {ai(k), a*(k)}. The field ip( h \ defined by 

(p( b \t, x, i) = U v i( b *j + h i) e~^ x+t - 1 ^ + i (b* - bj) e -%(^*+*™)] , (2.32) 

satisfies the equation of motion ( 12. ip and the boundary condition (12. 2p by construc- 
tion. The solution depends on a free parameter t m G M, whose physical meaning is 
clarified in the next section. 

In order to investigate locality, we compute the commutator 

[ipMfax^) , if {b) (t 2 ,x 2 ,i 2 )} = -2i £ e"^ 12 sinhirjj^U^U^ 

= -2i e "^ 12 smh(^ 12 ) Uiii U Ji\ > ( 2 - 33 ) 
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where X v = {i = 1, n : rji = 77}. Now, making use of the identity 

R ^ = b l^S k ~ iv) l us ^ u ~ x ] ix« = 2 E w ^ > ^ 2 - 34 ) 

one finds 

[ipVfaxuh) , ^ (6) (^,x 2 ,i 2 )] = -i e "^ 12 smhfe) R% 2 . (2.35) 

We stress that ( I2.35P is £ m -independent. Combining (I2.28P and (I2.35p . one gets for 
the commutator of the total field <p = + 

, (p(t 2 ,X2,i2)} = 

-~ [e(ti2 + ^12) + e(*i2 - ^12)] S hi2 - ~ [e(ii 2 + x i2 ) + e(t 12 - x 12 )] 5^(0) 

-ifl(f 12 - z 12 ) £ e ^-^)itg 2 + ^(-t 12 - x 12 ) £ e -^*-+^)i?g 2 (2.36) 

where V = V + UV-. Taking into account Xi j2 > 0, one has at space-like separations 
t\ 2 — X\2 < that \t\ 2 \ < |xi 2 | < |5i 2 |. These last inequalities imply the vanishing of 
(I2.36P at space-like distances, showing that ip is local. The canonical commutation 
relations (12.31 12 .4[) hold as well. We thus proved the following statement. 

Proposition 3. The boundary component ip^ b > defined above is such that the 
total field both canonical and local for any value of the parameter 

t m G BL 

Let us consider now more closely the time evolution of the system. The explicit 
form of the Hamiltonian H can be deduced from the identity 

[H , <p(t, x, i)} = -i(dt<p)(t, x, i) . (2.37) 

As expected, 

H = H {s) + H {b) , (2.38) 
where the scattering contribution has the standard form 

1 JL. r°° Ah 

H^ = -J2J ^\k\a*(k) ai (k). (2.39) 
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Concerning the BBS part H^ b \ eq. fl2.37p implies 



[ffW ,b j + b*}= vrjjibj + b*) , [H® , hi - &*] = -irjibj - b*) , j E X + . (2.40) 
which lead to the following expression 

H (b) = ~ E Vi {b) ~ bf) • (2.41) 

We see that is not at all the Hamiltonian of the conventional harmonic oscillator, 
but resembles that of the damped oscillator [32] . This fact has a relevant impact on 
the spectral properties of the total Hamiltonian and the time-evolution of the whole 
system, which are discussed in section 3. 

In conclusion, the main lesson from this section is that the BBS contribution ip® 
and the time evolution of the total field ip are fixed up to the parameter t m by the 
physical requirement of causality (local commutativity). 

2.4 Dual field, chiral fields and vertex operators 

This subsection is a collection of definitions of the basic structures needed for 
bosonization with BBS on T. The dual field (p, defined by 

d t <p(t,x,i) = -d x tp(t,x,i) , d x <p(t, x, i) = -d t (p(t, x, i) . (2.42) 

is given by 

<p(t, x, i) = p {s) (t, x, i) + ^ {b) (t, x, i) , (2.43) 

where 

ft s \t,x,i) = [°° dk£ ^ [oJ-C^e^l*'*-**) + OiC^e-^l , (2.44) 
J-oo 27Ta/2|A;| 

and 

p {b \t, x,i) = -^J2 U H [~ ( b *J + b j) Q~ V3{x+t ~ tm) + i (b* - bj) e-wf*-*^)] , 

(2-45) 

One easily verifies that <p> is a local field as well. We stress that ip and (p are however 
not relatively local, which is fundamental for bosonization. 
The chiral fields are given by 

Vi,R{t -x) = (p(t, x, i) + (p(t, x, i) , (p i>L (t + x) = ip(t, x, i) - (p(t, x, i) . (2.46) 
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The BBS contributions to (12.461) are 

if^t -x)=iV2j2 U*i { b ] - h) e-^(-*+*») , (2.47) 

+ x) = V2J2 % (b* + bj) e-^ +t ~^ , (2.48) 

and differently from the scattering parts do not oscillate in time, but vanish or 
diverge exponentially in the limit t — > ±oo. This is a first signal that the BBS 
give origin to a kind of instability in the theory, producing a complementary damp- 
ing/enhancement of (pf' R and ipf\ in time. We shall characterize this instability more 
precisely later on at the level of correlation functions. 

Once we have the chiral fields, we can introduce the family of vertex operators 
parametrized by £ = (a, r) £ M 2 and defined by 

V(t, x, i; C) = Zi Ki q(i; () : exp {iv^ [<J<Pi,n{t - x) + T(p i)L (t + x))} : , (2.49) 

where: 

(i) Zi £ R are normalization constants which are not essential in what follows; 

(ii) Ki are the so called Klein factors, which satisfy 

^i\^%2 ~l" ^12^11 = 2<5i 1 t a l , (2.50) 
and ensure [7] that the vertex operators (12.491) obey Fermi statistics provided that 

a 2 - T 2 = 2k + l, keZ; (2.51) 

(iii) the factor q(i; £) is given by 

q(i; C) = exp [iy/n (aQ itR - rQ ijL )] , (2.52) 

where 

1 f°° 

Qi,z = ~ d£ d^iAi) , Z = R,L. (2.53) 

^ J — oo 

are the chiral charges; 

(iv) finally, : • • ■ : denotes the normal product in the algebras A and B. 

The BBS contribution to the vertex operators (I2.49p is encoded in the chiral fields 
(J223 (MED- 



IO 



3 Correlation functions 



We start by emphasizing that the above analysis of the massless scalar field tp with 
off-critical boundary conditions on T shows that the theory admits two different 
regimes, which correspond to V+ = and V+ 7^ 0. Investigating the correlation 
functions, we will show now that these two regimes have very different physical 
properties. The considerations in the previous section are purely algebraic and do 
not depend on the representations of A and B. However, for the physical applications 
(correlation functions, transition amplitudes, etc.), we must choose a representation 
of each of these algebras. Up to unitary equivalence, the choice of the Fock repre- 
sentation T(B) of B is actually unique. This fact is essential for the case V+ ^ 0. 
In analogy with B, we take below the Fock representation J 7 (A) of A as well. The 
total Fock space is J 7 = J 7 (A) <8> ^(B), the vacuum being Q, = Q {s) <g> Q {b \ where 
fl^ and are the vacuum states in J 7 (A) and J-{B) respectively. For the vacuum 
expectation values in J 7 we adopt the following short notation 

(n,OQ) = (O) . (3.1) 

Let us discuss now the properties of the Hamiltonian (12.381) in this representation, 
focusing on the BBS contribution (I2.4ip . Since the *-operation in B is realized as 
Hermitian conjugation in J-(B), the BBS Hamiltonian H^ h > is a Hermitian operator. 
However, as shown in Appendix A, is not self-adjoint. In fact, its domain 
T>( b > C J-(B), where (12.401) hold, contains eigenstates with complex eigenvalues, 
implying that the time evolution of our system is not unitary. A direct proof of this 
statement is given in the next subsection. Analogous behavior, based on complex- 
valued Hamiltonians, has been used recently [33] in the description of dissipative 
quantum mechanics. 



3.1 The half-line 

We start our considerations with the half-line M + . Apart of being relatively simple, 
this case is fundamental because by means of (12. 12|) the derivation of the correlation 
functions on a generic graph V can be reduced to that on ]R + . 
On R + one has 

U = e 2ia , S(k) = , r) = X tan(a) (3.2) 

k — irj 

and the boundary condition (12. 2p takes the familiar form of mixed (Robin) condition 

(d x (p)(t,0) + v<p(t,0) = 0. (3.3) 
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S(k) has one BBS for rj > 0, the associated oscillator being generated by {b, b*}. 
For the basic two-point correlators one finds (£12 = £1 — £2, £12 = £1 + £2): 

W^i)M6))o = u( f i!; 12 ) + 2e( V )e T >&s- 2t ™\ (3.4) 
(<Pl(Zi)M&))o = u(nt 12 ) + 26( V )e-v^- 2t -\ (3.5) 

(MZi)<Pl&))o = 29(rj)[v+(-^ 12 ) - ie^ 12 ] + 2e(-rj)v^-^ 12 ) - «(Ki 2 ),(3.6) 
(Vl(6)M&))o = 2^( ?7 )K( ? 76 2 ) + ie-^ 12 ] + 2^-77)^(^2) - ^£12) , (3.7) 

with 

«(£) = -- + c) = — ln(|el) - Ut) , (3.8) 

7T 7T 2 

u±(0 = --e- f Ei(e±i6), (3.9) 

7T 

where e > and Ei is the exponential integral function. Finally, /j, is a infrared 
mass parameter (see e.g. [34J). The correlators (13.4113. 7p capture all fundamental 
properties of the theory and clearly show the difference between the two cases rj < 
and 77 > 0. 

For rj < the S'-matrix (13.21) has one antibound state. This case has been studied 
in detail in [34J. The correlators (l3.4H3.7p are invariant under time translations, the 
energy is conserved and <p has unitary time evolution. 

The presence of a BBS changes completely the situation for 77 > 0. The expo- 
nential terms in (13.4113. 7p . collecting in this case the contribution of the BBS on M.+, 
have two essential features: 

(i) The factor £ 12 in (13.41 13. 5p depends on ti + 1 2 , which shows that the theory is 
not invariant under time-translations. The origin of this symmetry breaking is the 
BBS contribution in the Hamiltonian (I2.39p . which does not annihilate the vacuum 

n, 

HQ = -^ v b* 2 Q ± 0. (3.10) 

Accordingly, for 77 > the energy of the system is not conserved, which signals a 
nontrivial energy flow crossing the boundary at x — 0. We stress that this is a 
purely boundary effect and that the energy momentum tensor in the bulk satisfies 
the continuity equation. 

(ii) For t m 7^ the correlators (I3.4ti3.7p are not invariant under the conventional 
time reversal transformation t — > —t. However, the combination of this operation 
with a specific time translation according to 

t^-t + 2t m , (3.11) 
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generates a symmetrj{§ of the correlators fl3.4H3.Tj) . It is implemented by 

T<p i)R (t - x)^ 1 = <fi>L (-t + 2t m + x), (3.12) 
Tip i>L (t + x)T- x = ip i>R (-t + 2t m -x), (3.13) 

T being an antiunitary operator which leaves invariant the vacuum, TQ = Q. 

A physical observable, which nicely illustrates the above features is the vacuum 
energy of the system. Let us derive first the vacuum energy density of the left and 
right movers separately. Like in the Casimir effect on the graph F [8J, this density 
is defined via point-splitting by subtracting the contribution on the line, namely 

6> z (0 = i lim [(tyz)(£i)(fy>s)tei) - ((^z)(6)(^z)(6))iine] , Z = L, R. 
2 £i,2->£ 

(3.14) 

Now, using (13. 4[ 13.5 j) one obtains 

£ L (t + x) = (6 L (t + x)) = V 2 e~ 2v{t+X ~ tm \ (3-15) 
£ R (t-x) = (e R (t-x)) = r } 2 e 2 ^- x -^. (3.16) 

For the total vacuum energy one gets 

POO 

E(t)= dx[£ L (t + x)+£ R (t-x)] = r)cosh[2r)(t-t m )]. (3.17) 
Jo 

As expected, the vacuum energy is time dependent: it decreases (dissipation) in the 
interval (— oo,i m ) and increases (enhancement) in (t m ,oo). E(t) has an absolute 
minimum at t = t m , which fixes the physical interpretation of the free parameter t m 
present in the general solution (I2.32[ I2.45[) . Obviously, t m indicates also the instant 
in which the vacuum energy flow trough the boundary changes its direction. Finally, 
being a reflection with respect to t m , the time reversal transformation (13. lip inverts 
the regimes of dissipation (t < t m ) and enhancement (t > t m ) of the vacuum energy. 

In order to understand better the above phenomena, it is instructive to study 
the time evolution of the vacuum state Q. Using the results of Appendix A, one can 
compute the expectation value 

(Q , e^fl) = (e itg ) = r =L= (3.18) 

which vanishes for t — > ±oo. In these limits therefore Q evolves in a state which is 
orthogonal to itself, which is essentially the result of [32] for the damped oscillator 



^Notice that the mapping (13.11[) is actually the reflection on R with respect to the point t„ 
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with 7] playing the role of friction. In this respect the vacuum energy densities (13. 16} 
13 . 1 5|) show the presence of positive friction in the left sector and negative friction in 
the right one. 

The existence of friction suggest also the violation of unitarity. In fact, deriving 
along the above lines the time evolution of the one-particle BBS state b*Q, one finds 

(e itH b*Q , e itH b*Q) = cosh(^) . (3.19) 

The time dependence of (13. 19[) implies the breakdown of unitarity. 

It is interesting at this point to investigate the time evolution of a generic coher- 
ent state 

Q c (z) =e zb *~* b n, z = x + iyeC, (3.20) 

generated by the boundary degrees of freedom. The main points of the computation 
are sketched in the Appendix A. The result is 

p — (x 2 +y 2 )[l+sech(rit)]— 2ixy ta,nh(qt) 

(n c (z) , e itH Q c (z)) = , (3.21) 

a/ cosher J 

showing that the vacuum Q and the coherent state tt c (z) behave in the same way 
in the limit t — > ±oo. 

The results of this subsection are collected in the following proposition. 

Proposition 4: The boundary bound state, which occurs for rj > in the mass- 
less local scalar field (p on the half-line R +J implies: 

(a) breaking of time translation invariance; 

(b) vacuum decay like in the case of the damped oscillator with friction rj; 

(c) nonunitary time evolution; 

In agreement with point (a), the dynamics of the system depends on the pa- 
rameter t m and exhibits two different regimes: the vacuum energy flow trough the 
boundary is outgoing for t < t m and incoming for t > t m . For simplicity we analyzed 
in this section the half-line, but we will show below that the above results hold for 
a generic star graph as well. 

3.2 Generic star graph 

The two-point correlation functions for a general star graph T are easily obtained 
combining the results on IR + with the diagonalized form f)2.12p of the S'-matrix. One 
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finds 



K*(£iW(6)>o = X) e ' (5a " 2tm) ^ + *i*i u (/'^) > (3.22) 
K,i#iW(6)>o = ^e-"^- 2 ^4t + ^ 2 ^a 2 ), (3.23) 

(^i,i?(6)^i 2 ,i(6))o = 

E K(-^i 2 ) - ie^ 12 ] 4t + E M-^i 2 )i?S 2 - <U , (3-24) 

(^i,L(6)^ 2 ,fi(6))o = 

E M^i 2 ) + ie-^ 12 ] itg 2 + E ^ + (^i 2 )4t - <U u(A*&0 • (3-25) 

In (13.221 13.23ft we recognize a t\ + ^-dependence. As expected, the presence of BBS 
(7-+ 7^ 0) implies the breakdown of time-translation invariance. Instead of (|3.18|) . 
one has on T 

r i r + 

(Q , e itH 0) = (e itH ) = == , (3.26) 

■\/ cosh (771) 

n + being the number of positive r]{. Equation (I3.26P shows the time evolution of 
the vacuum. For the vacuum energy density of the chiral fields on T, one gets from 
(I3T221 13T23D 



£%,L{t + x) = E v 2 ^ 



2 p -2ri(t+x-t rn )j^(v) 



£i,R,{t - x) = E v 2e 



2 p 2ri(t-x-t rn )j^(v) 



(3.27) 

which behave in much the same way as on the half-line. Finally, the total vacuum 
energy is 

Eit) = E VRti ] cosh[2r](t - t m )} . (3.28) 

In conclusion, the fields (p^z on F have the same features as their counterparts 
on the half-line. Proposition 4 is therefore valid on any star graph. 



4 Luttinger liquid with off-critical boundary con- 
ditions 

The quantum field theory on T, receiving much attention PQ-P2], is the Tomonaga- 
Luttinger (TL) model because it captures the universal features of a wide class of 
one-dimensional quantum many-body systems, called Luttinger liquids. 
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The dynamics of the TL model is defined by the Lagrangian density 



C = vpl{dt + vpd^x + iip* 2 {d t - v F d x )ip 2 - g+i^i^t + V^) 2 - g-Wx^x ~ ^l^f , 

(4.1) 



where ip a (t, x, i) with a = 1, 2 are complex fermion fields, vp is the Fermi velocity and 



showing that U(l) (g) U(l) is the (internal) symmetry group of the model. 

It is well known (see e. g. [351 US]) that the TL model on the line K. is exactly 
solvable by bosonization, the solution being expressed in terms of the massless fields 
ip and (p. On the graph T the situation is a bit more involved because of the 
additional interaction localized in the vertex V. At criticality, the solution via 
bosonization has been given in [71 E]. The basic physical observable, derived there, 
is the electromagnetic conductance matrix 



where G\- me is the conductance on the line and S is the scattering matrix which is 
constant because of scale invariance. Our main goal below is to derive the coun- 
terpart of (I4.4p away from criticality. For this purpose it is enough to consider the 
TL model in the special case when g + = —g~ = gn > and vf = 1. The classical 
equations of motion of this system, known as the Thirring model [37J, can be written 
in a simple matrix form 





(4.4) 



7 x d x )i/;(t : x,i) 



^g[ltJt{t, x, i) - j x J x {t, x, i))ip{t, x, i) , 



(4.5) 



where 




(4.6) 



and 



J u (t,x,i) — ijj(t,x,i)^ u ip(t,x,i) , 
is the conserved [/(l)-current. The f/(l)-counterpart is 



if) = ip*j t ■ 



(4.7) 



J u (t, x, i) = ip{t, x, i)ju75^(t, x, i) , 



75 



(4.8) 



which is also conserved. 



Another frequently used notation is #2 = .<?+ — g- and = g + + g_ . 
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4.1 Bosonization 

For solving the Thirring model, we set a > and 

1 1 

7p 1 (t,x,i) = -=V(t,x,i;(), ^ 2 (t,x,i) = -=V(t,x;C), (4.9) 
v /7T v 2n 



the vertex operators being defined by ( I2.49[) with ( = (a, r) and Q' = (r, a). In order 
to have canonical fermions we require 

<7 2 -T 2 = l, (4.10) 

implying a ^ ±r. We consider here the standard fermionic Luttinger liquid on T, 
but the results below have a straightforward extension to the anyonic case [13] as 
well. 

The quantum current J u is constructed by point-splitting according to 

J u (t, x, i) = - lim Z(e) \ip(t, x, i)^ v tp(t, x + e, i) + ifi(t, x + e, i)-y u ifj(t, x,i)] , 

(4.11) 

where Z(e) is some renormalization constant. The latter can be fixed [7] in such a 
way that 

J u (t,x,i) = -- — ■ — —=d v ip(t,x,i), (4.12) 
thus generating the [/(l)-phase transformations ( 14.21 ) 

[J t (t,a;,i) , ip(t,y,j)) = -S(x - y)5 ij ?p(t,xj) . (4.13) 
Because of (I4.12p the quantum equation of motion takes the form 

i{jtd t - Jxd x )ip(t, x, i) = - J 9 ^ 1 \ : (j t d t <P - l x d x <f) ip ■ (t, x, i) . (4.14) 

(a + r) 

Now, using the vertex realization ( 14.9}) of -0, one easily verifies that ( 14.141) is satisfied 
provided that 

r(a + r) = g. (4.15) 

Combining (I4.10p and (I4.15p . one determines a and r in terms of the coupling 
constant 

a = —== > , t = ■ , 4.16 

which fix the solution completely. 
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The quantum current J u is defined by point splitting as well. One has 

J u {t,x,i) = - lim Z(e) [ip(t, x, i)^ 5 i/;(t, x + e, i) + tp(t,x + e,i)^ 5 t{)(t,x,i)] 

(4.17) 

with 75 = —"itlx- A suitable choice [7] of the renormalization constant Z(e) leads to 

J u (t,x,i) = -- 1 d v <p(t,x,i) , (4.18) 

(j - r K/7T 



which is indeed conserved and generates the [/(l)-phase transformations (14.31) ac- 
cording to, 

r -5{x -y)S ij ip 1 (t,y,j) , a = l, 
[J t (t,x, l ),Mt,y,3)} = { S (x-y)5 ij Mt,y,j), a = 2. (419) 

It is worth mentioning that the boundary condition (I2.2p is transferred to a local 
condition on J u , even away from criticality. Indeed, (I2.2p implies 

n 

(dJ x )(t,0,i) + \J2 B iMt,0J)=0, (4.20) 

3=1 

where 

B = UBdU~ 1 , Bd = diag (tan(ai), tan(a 2 ), tan(a n )) . (4-21) 

A simple way of proving eq. (14.201) is to check it directly on the half-line (n = 1) and 
extend the result on T by means of (I2.10H2.131) . Using the duality relation (12.421) . 
the boundary condition (14.201) can be written also in terms of the [/(l)-current 

n 

(d x J t )(t, 0,i) + \J2 BijMt, 0, j)=0. (4.22) 

3=1 

We turn now to the important issue of symmetries. Clearly, the boundary condi- 
tions (12.21) influence the symmetry content of the model on V. The relation between 
symmetries and boundary conditions is implemented by the Kirchhoff 's rule, which 
must be imposed in the vertex of V on any conserved current in order to generate a 
time- independent charge from it. The matrix U in (I2.2p for instance parametrizes 
all boundary conditions which ensure the Kirchhoff rule for the energy-momentum 
tensor of tp and thus the time-independence of the relative Hamiltonian. Since dif- 
ferent conserved currents generate in general nonequivalent Kirchhoff's rules, one 
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may expect the presence of obstructions for lifting all symmetries on the line to 
symmetries on T. It may happen in fact that two Kirchhoff's rules are in contra- 
diction for generi^ boundary conditions. In this case one can preserve on T one of 
the corresponding symmetries, but not both of them. This is actually the case with 
the U(l) <S> [/(l)-group. 

Let us consider first the £7(l)-factor. Using (I2.14p . the associated current (I4.12p 
can be decomposed in scattering and BBS part, namely 

J u (t,x,i) = Jl s \t,x,i) + 4 b) (t,x,i) . (4.23) 

According to [7] 

n n n 

J24 s) Mi) = < > J^S^k) 1 < > ^V,-,- 1. Vj = l,..,n. 

i=l i=l i=l 

(4.24) 

For the BBS contribution one gets from (12.321) 

n n 

^jf(t,(U)=0 £^ = 0, VjGX + . (4.25) 

i=l i=l 

In Appendix B we show that the constraint (I4.24p on U actually implies the condition 
(I4.25P on U. Therefore the i7(l)-Kirchhoff rule holds if and only if the matrix U 
satisfies (14.241) . Analogously, for the [/(l)-Kirchhoff rule one gets 

n n 

Mt, o, i) = o u v = - 1 ' ( 4 - 26 ) 

i=\ i=\ 

Comparing the constraints (I4.24p and (I4.26P on the matrix U, we see that they 
cannot be satisfied simultaneously. Therefore, differently from the line, on T one 
must choose between the [/(l)-symmetry or its dual U(l). The corresponding sys- 
tems have different physical properties. Coupling minimally the model to an external 
electromagnetic field A u (t,x,i), one finds that (14.241) implies electric charge conser- 
vation. On the other hand, imposing (14.261) . the electric current J x does not satisfy 
the Kirchhoff rule and the electric charge is no longer conserved. However, equation 
(I4.26P implies by duality (I2.42p that the charge density satisfies 

n 

53 J t (t,0,i) =0, (4.27) 
i=i 

which is considered [11] to be the characteristic feature of a super- conducting junc- 
tion. Both (14.241) and (14.261) have therefore a precise physical interpretation. 

10 Excluding some exceptional boundary conditions in graphs with even number n = 2m of edges 
for which the system behaves as a bunch of m independent lines. 



19 



4.2 Conductance 

In order to derive the electric conductance tensor Gy, we couple the system to a 
classical external field A v (t,x,i) by means of the substitution 



d u + iA u (t,x,i) 



(4.28) 



in eq. (14. 5 p . The resulting Hamiltonian is time dependent and the conductance can 
be extracted from the linear term of the expansion of (J x {t, x, i))a v in terms of A u . 
This term can be computed by linear response theory [38] . Referring for the details 
to [7], one has 



n(l + 2g) 



(J x (t, x, i)) Av = {Jx(t, x, i)) + i / dT([H- mt (r) , J x (t, x, i)]} = 

J — oo 

n ft POO 

A x (t,x,i) / dr / d yAy{T,y,j){[d y ip{T,y,j) , d x (p{t,x,i)]) 

„•_ i J-oo JO 



(4.29) 

Let us consider now a uniform electric field"] E(t,i) = d t A x (t,i) which is switched 
on at t = to, i.e. A x (t, i) = for t < t Q . Using the Fourier transform of (I2.36p . one 
can compute the expectation value (14.29)) . The final result is 



(J x (t,0,i)) A „ 



n „ c 

Gline I 
3=1 J - ( 



did 

2^ 



where 



A x {u,j)e 



1 



— iu>t 



line 



AJu,i) 



dTe^A x (r,i) 



2tt(1 + 2p) ' 

The conductance can be extracted directly from (14.301) and reads 



(4.30) 



(4.31) 



Gij(u, t — to) 



G\\ 



line 



„ r>p( t ~ t o)'0 , \ 



T] + 100 



, t>t , (4.32) 



which is the off-critical generalization of ( 14.4)) . Notice that Gij(u, t — to) is in general 
complex, which leads to a nontrivial impedance [8]. The result (I4.32[) concerns the 
Thirring model, but the TL-conductance tensor is essentially the same, the only 
change being in the form of Gi ine which now depends [T5] on g± and vp. 
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We use the Weyl gauge A t = 0. 
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It is worth mentioning that (I4.32p satisfies the Kirchhoff rule 



Y,G ij {w,t-t o ) = 0, (4.33) 

3=1 

provided that the electric charge is conserved, i.e. (I4.24p holds. In fact, (I4.24p 
ensures the compensation of the first two terms in (I4.32p . As far as the Rij - 
contribution is concerned, combining (I4.24p with the definition (I2.26P one gets 

n 

53^ = 0, W V EV, (4.34) 

3=1 

which completes the argument. 

The off-critical conductance (14.321) keeps track of the two regimes V+ = and 
V+ 7^ of the theory. In fact, it is instructive to consider them separately. 

(i) V+ = 0: In this case the sum in the right hand side of (I4.32p runs only 
over the antibound states all of them producing damped oscillations in t — t . 
Therefore 

lim Gij (u, t - t ) = lim G^(u>, t - t ) = G hne [<% - SV,- (u)] , (4.35) 

t— >oo t — oo 

which gives the conductance, one will observe in this regime, long time after switch- 
ing on the external field. At a critical point S is constant and (14.351) reproduces 
(14.41) , which is claimed [U SJ EJ E] to describe some universal features of Luttinger 
junctions. 

(ii) V+ 7^ 0: Now the sum in the right hand side of (14.321) involves terms with 
7] > 0, which give origin to oscillations whose amplitude is growing exponentially 
with t — 1 . 

It is worth stressing that presently we are not aware if the response of real-life 
quantum wire junctions to an external electric field can be of the type (i) and/or 
(ii). It is nice however that the predicted behavior of the conductance in these two 
cases is quantitatively different, thus providing a clear experimental signature. 

The result (I4.32p can be applied also at the level of effective theory. An example, 
which frequently appears in the literature jU EJ E], is the Y-junction with the geom- 
etry displayed in Fig. 2. One has three external half lines Ei and a ring composed 
of three compact internal edges Ii. A magnetic flux (p is crossing the ring. The 
complete field theory analysis of the Luttinger liquid on a graph with this geometry 
is very complicated problem. One approximate way to face this problem could be 
to use the star product approach [UJ or the "gluing" technique (121 S3] for deriving 
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Figure 2: A graph with 3 external and 3 internal edges. 



the 3x3 scattering matrix relative to the external edges. This ^-matrix can be used 
for developing a simplified model with one effective off-critical junction for which 
(pQ2D applies. 

In conclusion, we emphasize that the results of this subsection are obtained in 
an abstract setting, where we attempt to describe the physics of a quantum wire 
junction by a Luttinger liquid with specific off-critical boundary conditions at the 
vertex of a star graph, approximating the junction. Further investigations are needed 
for clarifying both the plausibility of these assumptions and the applicability of the 
result (I4.32p to realistic quantum wire junctions. 

5 Outlook and conclusions 

We started this paper by developing the theory of a massless scalar field ip with 
off-critical boundary conditions on a star graph. The vertex of the graph is repre- 
sented as a point-like defect characterized by a scattering matrix S(k). In general 
S(k) admits bound and antibound states. This fact determines two physically dif- 
ferent regimes of the theory. In absence of bound states the field <p has unitary 
time evolution and the energy is conserved. The behavior of the system changes 
completely if bound states are present. Each such a state generates an oscillator 
in the spectrum of ip, whose contribution is fixed by local commutativity up to a 
common free parameter t m . These oscillators break the invariance under time trans- 
lations and drive the system out of equilibrium. Accordingly, the vacuum energy is 
time dependent: it decays exponentially for t < t m and grows at the same rate for 
t > t m . The two regimes are related by time reversal and correspond to a nontrivial 
outgoing and incoming vacuum energy flows. Recalling that the traditional way [39] 
for driving a system out of equilibrium is to couple it with a "bath" of external 
oscillators, we discovered above that such oscillators are automatically present for 
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S(k) with bound states. Such boundary conditions provide therefore an intrinsic 
mechanism for constructing non-equilibrium quantum systems. It is worth stressing 
that this mechanism is based on purely boundary effects and is completely fixed by 
the fundamental physical requirement of causality. 

It turns out that the behavior of the Luttinger liquid on T is affected by the 
above phenomena. This is true in particular for the electromagnetic conductance 
G. Deriving an explicit formula for G, we have shown that it depends on the time 
elapsed after switching on the external electric field and develops oscillations with 
both exponentially growing and decaying amplitudes, corresponding to the bound 
and antibound states respectively. 

Being interested in the Luttinger liquid, we investigated the massless scalar field 
on T, but the BBS phenomena described in this paper take place also in the massive 
case m > 0, provided that the boundary parameter 77 satisfies 77 > m. In the range 
< 77 < m the BBS do not produce instabilities [40J . 

Let us comment finally on some possible generalizations of the above results. It 
will be useful for the applications to extend our framework to boundary conditions 
which violate time reversal [21 IS]- This more technical then conceptual problem has 
been addressed recently within the above framework in [T5]. The study of the rich 
spectrum [T3] of effects away from equilibrium is essential for the experiment. More 
general systems of wires with several junctions attract also much attention [4"I]-|45]. 
Networks involving compact edges, like that in Fig. 2, are expected to have richer 
BBS content and represent a challenging open problem. 

We stress that the phenomena, investigated in this paper, illustrate actually very 
general features of systems with boundaries. In fact, the BBS mechanism, described 
above, has a straightforward extension |46] from point like to higher dimensional 
boundaries and admits therefore other potential applications. 
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A BBS Hamiltonian and vacuum decay 

We demonstrate first that the BBS Hamiltonian f)2.40p is Hermitian but not self- 
adjoint. It is enough to consider the case of one BBS, i.e. 

H® = l -r]{b 2 -b* 2 ) , 77 >0. (A.l) 
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This Hamiltonian generates the correct time evolution of ip( b \ provided that and 
b ± b* have a common and invariant domain D' 4 ' where the commutation relations 

[H^ ,b + b*}= irj(b + b*) , [H® ,b-b*] = -irj(b - b*) , (A.2) 

hold. We show now that is Hermitian but not self-adjoint on P 6 '. In fact, let 
$ e £> (b) be an eigenstate of H^> with eigenvalue E, i.e. = E&. Suppose 

that BgI, otherwise there is nothing to prove. Using flA.20 one easily verifies that 
is an eigenvector of H^> with eigenvalue E + i?7, which is complex 
since rj > 0. This concludes the argument. 

Our next step is to sketch the derivation of equations (I3.18P and (|3.2ip . We start 
with f)3.18p . observing that 

(Q,e itH Q) = (Q,e^ 2 - b ^Q), (A.3) 

where 

[b , b*} = 1 . (A.4) 

Setting M = —2 — Ab*b one immediately verifies that 

[b* 2 , b 2 } = M , [J\f, b* 2 } = -8b* 2 , [M , b 2 } = 8b 2 , (A. 5) 

showing that {b 2 , b* 2 , A/"} constitute a generator basis for a realization of the SU(1, 1) 
Lie algebra. Therefore, one can rewrite the SU (1,1) group element e a ^ b * ~ b ) as a 
product of exponentials of the algebra generators, i.e. 

e *(b* 2 -b 2 ) = e f(a)b* 2 e g(a)b* e h(a)M ^ ( A _g) 

where /, g, h are some functions of the real variable a, satisfying 

f(p) = g(P) = h(0) = 0. (A.7) 

We are particularly interested in h, because 

(fi,e¥^ 2 - 62 ^)=e~ 2 Kf). (A.8) 

From group theory we know that f,g,h satisfy first order coupled differential equa- 
tions. They are easily obtained by deriving both sides of eq. ( 1A.6I) with respect 
to a and identifying the coefficients of b* 2 , b 2 and M after shifting to the right the 
remaining exponential part. In this way one gets 

df 



da 
dfl- 

da 
dh 

da 



1-4/ 2 , 
8/<7 - 1 , 

/• (A.9) 
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The solution of this system with initial condition ( \A.7\i is 



f(a) = - tanh(2a) , g(a) = — -sinh(4a) , h(a) = - In [cosh(2a)] , (A. 10) 

which, combined with ( 1 A . 8 1) . implies ( I3.18p . 

In order to derive (I3.2ip one can use the identity 

(e zb *-* b fl , e a(b * 2 ~ b2) e zb *-* b Q) = e^ 2 (fi , e^e /(a)fe2 e 9{a)b * 2 e hia)M e zb * Q) , (A.ll) 

which follows from (IA.4j) . flA. 6|) and bQ = 0. The strategy for computing the right 
hand side of flA.lip is simply to move the exponents involving M and b 2 and b to 
the right by using the commutation relations among the operators in the exponents. 
After some algebra one obtains (I3.2ip . 



B BBS contribution and Kirchhoff 's rules 

We show here that 

n n 

J2U tJ = ±l $^ = 0, (B.l) 

for all j such that a, ^ 0, ±~ in (12.111) . Rewriting 

U d = U- l UU (B.2) 



in the form 

and using that 
one gets 



n 



k=l k=l 

{U d ) kj = 6 k ,e 2ia " (B.4) 



k=l 

Summing over % from 1 to n one obtains from (1B.5P 

n n n 



i=l i=l k=l 
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Combining the hypothesis in (IB. 1|) with ( IB.6P and using U l = U, one finds 

n n 

]T^e 2 ^=±5> fcj , (B.7) 

i=l k=l 

or, equivalently 

n 

(e 2 ^Tl)^% = 0, (B.8) 

i=l 

which proves the statement (IB. II) . 
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